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We consider trial wavefunctions exhibiting SU(A^) symmetry which may be well-suited to grasp the physics 
of the fractional quantum Hall effect with internal degrees of freedom. Systems of relevance may be either 
spin-unpolarized states (K = 2), semiconductors bilayers (K = 2,4) or graphene (K = 4). We find that some 
introduced states are unstable, undergoing phase separation or phase transition. This allows us to strongly 
reduce the set of candidate wavefunctions eligible for a particular filling factor. The stability criteria are obtained 
with the help of Laughlin's plasma analogy, which we systematically generalize to the multicomponent SU(A") 
case. The validity of these criteria are corroborated by exact-diagonalization studies, for SU(2) and SU(4). 
Furthermore, we study the pair-correlation functions of the ground state and elementary charged excitations 
within the multicomponent plasma picture. 
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I. INTRODUCTION 

Soon after the discovery of the fractional quantum Hall ef- 
fect (FQHE), 1 Laughlin successfully described the underlying 
strongly-correlated electron liquid with the help of a simple 
trial wavefunction. 2 The reasons for the success of this ap- 
proach were twofold: first, the calculated energy of this state 
is lower than that of charge-density waves or Wigner crystals, 3 
which are natural candidates for the ground state within the 
partially filled lowest Landau level (LL) due to the quenched 
kinetic energy. 2 The second reason for its success is the fact 
that Laughlin's wavefunction is in very sharp agreement with 
exact-diagonalization studies^ 

A powerful tool in the understanding of Laughlin's wave- 
function is the quantum-classical analogy, in which its prob- 
ability is interpreted as the (classical) Boltzmann weight of 
a two-dimensional (2D) one-component plasma (2DOCP). 2 
Most strikingly, this analogy shows that Laughlin's trial wave- 
functions have no free parameter to be optimized by any vari- 
ational calculation. 

Laughlin's original proposal was concerned with only one 
single species of fermions, namely spin-polarized electrons. 
In spite of its success, this is at first sight a very crude as- 
sumption in view of the relatively weak (effective) Zeeman 
effect when compared to the leading energy scale set by the 
Coulomb interaction e 2 /e/g, in terms of the magnetic length 
Ib = y/ft/eB. Indeed the latter is almost two orders of mag- 
nitude larger than the bare spin splitting for typical mag- 
netic fields of B ~ 10 T. In order to account for an internal 
SU(2) spin symmetry, Halperin proposed a generalized trial 
wavefunction, 5 which includes Laughlin's as a special case. 
The latter may indeed be viewed as a Halperin wavefunction 
with a spontaneous ferromagnetic spin ordering. 6 

Halperin's SU(2) wavefunctions have been a first step in 
the understanding of general multicomponent systems. In the 
case of bilayer quantum Hall systems, the same wavefunctions 
may be applied if one supposes a complete polarization of the 
physical spin and if one interprets the two layer indices as the 
two possible orientations of a pseudospin. 6 - 7 However, the hy- 
pothesis of complete spin polarization is a priori as feably jus- 



tified in bilayer as in monolayer quantum Hall systems, again 
due to a relatively weak Zeeman effect. A more appropriate 
approach is therefore one that takes into account the internal 
SU(4) spin-pseudospin symmetry. Such approaches have in- 
deed been proposed in the description of ferromagnetic states 
when the LL filling factor v = n e //«B, in terms of the elec- 
tronic, n e u and the flux, ris = eB/h, densities, respectively, is 
1,2, or 3& 

Another example of a multicomponent quantum Hall sys- 
tem is graphene, where the internal SU(4) symmetry is 
due to the physical spin accompanied by a twofold valley 
degeneracy . 10 i n i 12 In contrast to the abovementioned bilayer 
quantum Hall systems, where the pseudospin symmetry is 
explicitely broken because of the difference between intra- 
and interlayer Coulomb interactions, the SU(4) symmetry is 
almost perfectly preserved in graphene from an interaction 
point of view - a possible (valley) symmetry breaking may 
be due to lattice effects, which are suppressed by the small 
parameter a/ls, where a = 0.14 nm is the distance between 
nearest-neighbor carbon atoms in graphene, as compared to 
Ib = 26 VB[T] nm . 13 ' 14 ' 15 i 16 ' 17 In order to describe a possible, 
yet unobserved, FQHE in graphene, taking into account the 
appropriate form of the interaction potential, 1318 exact diag- 
onalization studies have been performed in the framework of 
an internal SU(2) valley symmetry, 19 - 20 as well as in a SU(4) 
composite-fermion appro achM^ 

More recently, two of us have proposed a generalization 
of Halperin's wavefunctions to K components, i.e. systems 
with an internal SU(^T) symmetry, in order to describe a pos- 
sible FQHE in ^-component systems, namely graphene with 
K = 4. 23 Here, we investigate the stability of these wavefunc- 
tions from two complementary perspectives - first, we derive 
stability criteria within a generalized plasma picture. This 
analogy allows one to interpret the SU(-ST) Halperin wavefunc- 
tions in terms of K correlated 2DOCP and to describe in a 
compact manner their ground-state properties as well as the el- 
ementary excitations with fractional charge. In a second step, 
we corroborate the validity of the generalized plasma picture 
with the help of exact-diagonalization studies. 

After a brief review of Laughlin's plasma analogy (Sec. II), 
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we generalize the plasma picture to ^-component systems in 
Sec. III. In Sec. IV, we derive general stability criteria within 
the plasma analogy, on the basis of which we discuss the sta- 
bility of specific SU(2) and SU(4) wavefunctions. We com- 
plete this paper with a discussion of ground-state properties, 
such as sum rules for the pair correlation functions (Sec. V), 
and fractional charges of quasiparticle/-hole excitations (Sec. 
VI). 



It is evident from the last equation that the variational param- 
eter must be positive - otherwise one would have to treat with 
unphysical negative densities. Notice that from the wavefunc- 
tion point of view, m < is not physical because it violates the 
analyticity condition for wavefunctions in the lowest LL. This 
point, which may seem obvious, is worth being emphasized 
and being recalled in the following sections when the plasma 
picture is generalized to more components. 



II. LAUGHLIN'S PLASMA ANALOGY 



In order to describe a correlated electron liquid to ac- 
count for the FQHE, Laughlin proposed the TV-particle trial 
wavefunction 2 - 
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where Zk — Xk + iy k denotes the position of the k-th electron in 
the complex plane. Here and in the following, we set the mag- 
netic length ij e 1, for notational convenience. The form of 
this trial wavefunction is solely determined by the analyticity 
condition for the lowest LL - i.e. all single-particle states are 
of the form z c exp(-z 2 /4), where I is a positive integer - and 
by symmetry considerations. In order to have a translational 
and rotational invariant state and thus an incompressible state 
with no gapless Goldstone mode, the wavefunction may only 
depend on the relative distance Zk - zi of the k-th and the /-th 
particle. Furthermore, fermion statistics for electrons requires 
the exponent m to be an odd integer, which is the only varia- 
tional parameter of Laughlin's wavefunction (Q]l. 

However, the parameter m turns out to be fixed by the elec- 
ton density, or else the filling factor, v = 1/m, as Laughlin 
showed with the help of a plasma analogy. 2 Indeed, one may 
interpret the modulus square of the wave function as the Boltz- 
mann weight 
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of a classical system, namely a 2DOCP described by the clas- 
sical Hamiltonia n 2 ' 24,25 



<H N = -mJ]ln\z k -zi\ + J]^f, 



(3) 
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where one has set somewhat arbitrarily the inverse "temper- 
ature" /? E= 2. Notice that the true temperature does not in- 
tervene in the analysis because the system is placed at T — 0. 
The first term describes 2D interacting particles of charge y[m, 
whereas the second term may be interpreted as a homoge- 
neous background of charge - 1 / (jellium). The minimiza- 
tion of the classical Hamiltonian corresponds, via the relation 
©, to a maximal quantum probability of the original quantum 
system of electrons within the lowest LL. The classical ground 
state of the Hamiltonian (|3), however, is obtained when the 
plasma particles of charge -ym are fully neutralized by the 
background, i.e. when 
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III. PLASMA PICTURE FOR SU(K) HALPERIN 
WAVEFUNCTIONS 



Based on Halperin's idea to write down a Laughlin-type 
wavefunction for a two-component quantum Hall system, in 
order to take into account the spin degree of freedom, 5 two of 
us have proposed a SU(^T) generalization for a /^-component 

system,— 



■SU(K) 
m- i ,...,m K ;n ij 



K Ni 



x 0"" c/ ' x exp 

ft >i a 



1=1 ki<k 
K N, Nj 

iiiiii. . ■ 

i< j kj = l kj=l 



(5) 



There are K different types of electrons [denoted with super- 
script (/)] with inter- (« ;j ) and intra-component (m,) quantum 
correlations, and zl is the complex position of £,-th electron 
of type i = 1, K. The lowest-LL analyticity condition im- 
poses that all exponents, m, and «y, must be integers. Fur- 
thermore, nti must be odd in the case of fermions. Apart from 
K = 2, discussed by Halperin, 5 K = 4 wavefunctions may 
be physically significant in the case of bilayer quantum Hall 
systems and graphene. In the former example, the internal de- 
grees of freedom do not only contain the physical SU(2) spin 
(T, I), but also a layer index, which may be mimicked by an 
additional SU(2) isospin (+, -). There are thus four internal 
states, 1 = (T, +), 2 = (T, -), 3 = (I, +), and 4 = Q, -). In 
the case of graphene, an isospin (+, -) must be introduced in 
order to account for the two-fold valley degeneracy. Wave- 
functions similar to those in Eq. (O have been proposed by 
Qiu et aJ. for multilayer quantum Hall systems, 26 by Morf as 
potential candidates for the FQHE hierarchy states, 27 and by 
Yang et aJ. in the study of a possible FQHE in graphene.— 

Again, the starting point of the plasma analogy is Eq. (f2]), 
and one associates the new Hamiltonian "K^ with a physi- 
cal system. In the case K — 2, this system has been inter- 
preted as a generalized plasma, which consists of K different 
particle types (each of which corresponds to a different elec- 
tron type in the original quantum system) plus a neutralizing 
background. 6,28 With the identification (O, one obtains for the 
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wavefunctions (0 the classical Hamiltonian 
<H N = -t^t^k-^-t^Z^-^' 

i=l ki<l, i<j ki,kj 

K N, l z ®|2 
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(6) 



Here the first term represents a sum over K 2D interaction 
terms for (/)-type particles of charge yfmi, whereas the sec- 
ond one takes into account interactions between particles of 
different type, (/) and (J). However, this generalized plasma 
does not satisfy the charge superposition principle 29 unless 
«y = -Jmimj, which is a rather special case. 

Instead of one single plasma of K types of particles, it 
seems therefore more appropriate to interpret this generalized 
plasma in terms of K different 2DOCPs (one for each type 
of electrons) with correlations between them. For this pur- 
pose, we introduce the continuum limit, in which the density 
for particles of type (/) (electrons or plasmatic particles) is 
Pi( r ) = Yjki £( r ~~ r k,)- ln order to distinguish the resulting 
Hamiltonian from the original discrete one, we supress the N 
subscript in Eq. (0, and one obtains the energy functional 
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Here, Mk is the symmetric exponent matrix, with «,y = riji 
and Ha = mi 23 and £2 is the surface occupied by the plasma. 
Similarly to the one-component case, the configurations with 
maximal probability [Eq. (fJJ] are obtained by minimizing <H 
with respect to all densities. The stationary points are found 
at 
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In order to have a minimum, the Hessian matrix 
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which is identical to the exponent matrix M K up to a positive 
constant, needs to be positive, i.e. have positive eigenvalues. 
One may interpret Eq. ([8]l as the stationary point of a 2DOCP 
of (0-type particles, whereas the positions of all other parti- 
cles of type (j + i) are fixed and constitute a quasi-static im- 
purity potential felt by the (/)-type particles. For the 2DOCP 
of this type, the interactions between all other types of parti- 
cles yields only an unimportant constant, with respect to the 
pi derivative. In this sense, one may indeed interpret the sys- 
tem as K correlated 2DOCPs rather than a single plasma of K 
different types of particles. 



In the same manner as for a single 2DOCP, Eq. © is sat- 
isfied when each of the K plasmas exhibits quasi-neutrality, 25 
but now contributions from the impurities have to be taken 
into account, 
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Here, l/2n on the r.h.s. of Eq. (1 1 Oat represents the neutral- 
izing background, as for a single 2DOCP [Eq. (0J], The sec- 
ond term in Eq. (llOat represents the contributions from type- 
U + particles due to inter-component correlations. One 
notices that Eq. (llObb is the matrix generalization of Eq. ©. 
This result was previously derived by counting the zeros of 
the wavefunction (|5]l. 23 Invertible matrices yield a unique so- 
lution with all densities being uniform, p,(r) = p,, as it is the 
case for U(l) Laughlin's liquid. The case of non-invertible 
matrices will be discussed in the next section. 

Unlike the U(l) case, fixing the total filling factor vj = 
v\ + ... + vk does not determine uniquely the exponent ma- 
trix. It has been pointed out in Ref. 23 that several candidate 
wavefunctions may give rise to a FQHE at the same filling 
factor vt, especially in the case of larger internal symmetry 
groups. Moreover, even if one fixes all component filling fac- 
tors Vj = 2npi, the wavefunction is not unambiguous. As 
an example, we consider the SU(2) Halperin wavefunctions 
{m\,m.2,n), with (i) m\ = 1112 = 3, n = 1 [(331) wavefunc- 
tion], and (ii) nt\ = nix = 1, n = 3 [(113) wavefunction]. Both 
wavefunctions describe a situation with vi = Vz = 1/4 and 
have been considered in the past within the study of a possible 
unpolarized vj = 1/2 stated Indeed, an even-denominator 
quantum Hall state has been observed at v = 5/2 and v = 7/2, 
in the first excited LL. 31 However, it is strongly unlikely that 
this state is spin-unpolarized 32 ' 33 i 34 and more sophisticated 
theories, in terms of a Pfaffian state, need to be invoked to 
account for a spin-polarized stated As we will show below, 
stability conditions related to Eq. © allow one, in the case of 
unpolarized states, to discriminate between (mmri) and (nnm) 
wavefunctions, for m > n. 



IV. STABILITY 

In order to obtain a stable state, the stationary point ob- 
tained from Eq. dHJ must be a minimum, i.e. the Hessian 
matrix oc Mk in Eq. © must be positive. Otherwise, the plas- 
mas would not be stable, or else in the state of lowest energy. 
Hence, a first stability condition imposes that all eigenvalues 
A, of Mk must be positive. 

There is indeed a limiting case for which one, or more, 
eigenvalue(s) is (are) zero. Because the potential is quadratic, 
the minimum point becomes now a line of minima, which cor- 
respond to different ground states. If at least one eigenvalue is 
zero, the exponent matrix Mk is no longer invertible, and all 
densities may not be fully determined from Eq. (llObb . 
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This point may alternatively be interpreted in terms of 
SU(X) ferromagnets - different states of equal energy may, 
e.g., occur at various combinations of two (or more) fill- 
ing factors v; and Vj, although the sum Vy = v, + v, is 
fixed. In this case one may introduce a pseudo-spin opera- 
tor S*j = N{vi - Vj)/2, which can possibly take all values in 
between -N/2 < S|. < N/2. The simplest example of such 
a case is the Laughlin wavefunction with an internal spin de- 
gree of freedom [a Halperin (mmm) wavefunction], which is, 
for odd m, completely antisymmetric in its orbital part. For 
fermions, the spin wavefunction must therefore be completely 
symmetric and represents thus a SU(2) ferromagnet. If the to- 
tal spin is oriented along the positive z direction, all electrons 
reside in the upper spin branch (vi = 1 /m, V2 = 0). In the ab- 
sence of a Zeeman effect, this state has the same energy as the 
one with a total spin in the — z direction (vi = 0, V2 = 1/m), 
as well as any intermediate state with v\ + V2 = 1/m. In the 
general SU(ZQ case, the ferromagnetic properties are deter- 
mined by the rank r of the matrix M K . Indeed, if r < K and 
one introduces common fixed filling factors Vu for the relevant 
components, one may describe the resulting state by a SU(r) 
Halperin wavefunction with an invertible exponent matrix M r 
with additional pseudo-spin degrees of freedom for the com- 
ponents the density of which remain undetermined^ 

As mentioned earlier, some exponent matrices can lead to 
negative density solutions for Eq. ( llObl i. A second class of 
stability conditions needs to be imposed in order to prevent 
this unphysical situation, which may occur even in the case of 
a positive matrix Ma- 
in order to illustrate the two conditions, we discuss some 
specific examples for different K. The case K — 1 has already 
been presented above. 



A. The case K = 2 

We first study Halperin's wavefunction (ni\,mi,ri) for the 
SU(2) case, which is described by the exponent matrix 



M K=2 = 



nt\ n 
n ni2 



Even if all exponents are positive, as required by the lowest- 
LL analyticity condition, the eigenvalues A ± and the filling 
factors V1/2 are not necessarily so, 



m\ + ni2 + V( m i — 1112) 2 + 4n 2 



V\ = 



ni2 — n 
m\ni2 — n 2 



2 

V2 = 



nt\ — n 
m\ni2 — n 2 



(Ha) 
(lib) 



In order to obtain only positive eigenvalues (first stability con- 
dition), one needs to require 

m\m2 — n 2 > 0. (12) 

The case m\ni2 = n 2 corresponds to a situation of a non- 
invertible matrix of rank r — 1 . Because of Eq. ( fT2l . positive 
densities (filling factors) are found from Eq. ( 1 1 lbb only for 
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FIG. 1: Stability of the (73fi) wavefunction. Both filling factors, vi (long 
dashed line) and V2 (short dashed line), and the A- eigenvalue (solid line) of 
the (73«) wavefunction are plotted as a function of n. For n < 3 (part I) all 
quantities are positive and the corresponding state is stable. In part II, one of 
the filling factors is negative. Part III exhibits positive filling factors but the 
state is still unstable because of the negative eigenvalue; the system eventually 
undergoes a phase separation. 



(a) 



(b) 



mi > n and 



m,2 > n, 



(13) 



FIG. 2: Sktech of the plasma corresponding to the (73m) wavefunction, 
for (a) n < 3 and (b) 1 < n. (a) Intra-component are stronger than inter- 
component repulsions (n < m\,mz), Type-(l) particles (black) are there- 
fore on the average surrounded by type-(2) particles (grey). This yields a 
stable state of two homogeneous interpenetrating plasmas, (b) If the inter- 
component is stronger than the intra-component repulsion, the plasmas have 
a tendency to phase-separate to minimize the number of neighbors from dif- 
ferent types. 



which we thus need to impose as a second class of stability 
conditions. One furthermore notices form Eqs. ( TTZb and ( fT3b 
that the only states with a non-invertible exponent matrix of 
rank r — 1 are the ferromagnetic Laughlin states (mmm) dis- 
cussed above. 

The final stability criterion Eq. ( TT3l ) for SU(2) wave- 
functions has a compelling physical interpretation: intra- 
component must always be strongerthan inter-component cor- 
relations. Within the plasma picture, this may also be under- 
stood from Figs. [T][2{a), and|2|b). For illustration, we con- 
sider the (73n) wavefunction, where n is left as a variable, 
which we treat in a rather artificial manner as a continuous 
variable in the following discussion. It is evident that only in- 
teger values may be taken into account for physical candidate 
wavefunctions. 

Fig. Q]shows the plot of both component filling factors and 
the lower eigenvalue, /!_. This graph can be split into three 
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distinct parts 1,11, and III. 

• Part I (n < 3). Each type-(l) particle carries a charge 
V7 and is affected by those from the type-(2) plasma 
through a charge coupling of n/ V7 < V7, which may 
be interpreted as constituting a quasi-static impurity 
distribution interacting with type-(l) particles. Alter- 
natively, one may concentrate on type-(2) particles with 
charge V3, which see type-(l) particles as a distribu- 
tion of charge nj V3 < V3 impurities. Therefore both 
types of particles are more strongly repelled by those 
of their own species than by particles of different type. 
One thus obtains a stable homogeneous mixture of two 
plasmas, which is shown in Fig. [2ja). 

• Part II (3 < n < 1). Although type-(l) particles are still 
more strongly repelled by those of their own species 
than by type-(2) particles (n/ V7 < V7), this is not the 
case for type-(2) particles. Because their inter-species 
repulsion is now weaker than that which they experi- 
ence from type-(l) particles (n/ V3 > V3), they pre- 
fer to gather rather than to be mixed to type-(l) parti- 
cles. This indicates a tendency to phase-separate, and 
the effect manifests itself in an unphysical negative fill- 
ing factor v\. The divergence of the filling factors at 
the artificial value of n c = V21 ^ 4.6 in Fig. Q~|is due 
to the vanishing of the eigenvalue Above n c , the 
original minimum of the energy functional (0 evolves 
into a saddle point, and the first stability condition of 
non-negative eigenvalues is no longer satisfied. Indeed, 
one notices that the filling factors interchange their role 
- although the repulsion between type-(l) particles is 
stronger than that between type-(2) particles and one 
would therefore intuitively expect that v\ < vi, one 
finds vi > V2 for n > n c . 

• Part III (7 < n ). Above n — 7, the inter-component 
repulsion is stronger than that between particles of like 
type. The phase separation between the two plasmas, 
which we have alluded to in the discussion of part II, 
is well-pronounced. Due to this strong inter-component 
repulsion, the interface between the two plasmas needs 
to be minimized, and this results in an inhomogeneous 
state of two spatially separated plasmas, as shown in 
Fig. HJb). Should the partial densities not be fixed, i.e. 
particles could flip from 1 to 2 state, it is clear that the 
system would favor a distribution where only one type 
of particles would remain. The interface between the 
two plasmas disappears. 

Fig. [3]shows the stability graph for the a (mmri) wavefunc- 
tion (here with n = 3). In this case, the eigenvalues ( 11 1 al l and 
the component filling factors ( II lbl i become 



m ±n 



and 



Vi = v 2 



m + n 



(14) 



respectively. Although the component filling factors remain 
positive for all choices of n, the eigenvalue A_ becomes neg- 
ative for n > m, where one would expect a phase separation 




FIG. 3: Stability of the (33«) wavefunction. Equal filling factors, vi and V2 
(long dashed line), and the A- eigenvalue (solid line) of the (33;i) wavefunc- 
tion are plotted as a function of n. For n < 3 (part I) the state is stable, as in 
Fig. \T\ There is no corresponding part II, since the densities never vanish and 
in part III the plasmas tend to phase separate. 



between the two plasmas, as in the case of the (73n) wave- 
function [Fig. a)]. The critical value n = m corresponds to 
the Laughlin case with a SU(2) ferromagnetic spin wavefunc- 
tion, discussed above. For the case n < m, both trial wave- 
functions, (mmri) and (nrirn), are valuable candidates for the 
description of a potential FQHE at vj — 2/(m + n) if only the 
symmetry considerations for trial wavefunctions in the lowest 
LL are taken into account. However, the plasma analogy indi- 
cates clearly that only one of the two wavefunctions, namely 
(mmri), yields a stable physical state. At half-filling, e.g., the 
only SU(2) Halperin wavefunction which might yield a stable 
FQHE state is (331), whereas (1 13) corresponds to an unstable 
plasma, which is not evident from wavefunction calculations 
alone. 30 



B. The case K = 4 

We now consider generalized Halperin wavefunctions with 
an internal SU(4) symmetry. We restrict our studies to a par- 
ticular subset of the latter, noted (m\m2m\m2,n e n+nJ), the 
corresponding exponent matrices of which may be written 

„„23 



'm\ n e n + n e 
n e mi n e «_ 
n+ n e m\ n e 
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mi 



(15) 



If applied to graphene, those correlation coefficients imply 
that one treats all intervalley components (n e ) on the same 
footing and intravalley ones separately («+,«- for + and - val- 
ley, respectively). This is an even more natural assumption in 
the case of bilayer quantum Hall systems in semiconductor 
heterostructures, where interlayer correlations (described by 
the exponents n e ) are weaker than intralayer ones (n+ and n_ 
which couple the different spin orientations within the + and 
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- layer, respectively). Moreover, some intracomponent corre- 
lations are fixed to the same value for explicit calculation of 
the eigenvalues A and filling factors v to be carried out. How- 
ever, we will only settle here the conditions for all quantities 
to be positive (first and second stability arguments), and these 
conditions are satisfied if 

' m\ > n + 
m 2 >n- 
ni2 + «- > ln e 
mi+n + > 2n e . 

The case where one of the two first inequalities, or the two 
last, turn into an equality corresponds to matrices of rank 
r < 4. For r — 1 Eq. (TT&t becomes an equation set and, once 
again, the only stable state is of Laughlin-type (mmmm, mmm) 
with SU(4)-ferromagnetic ordering. The stability criteria [Eq. 
([ToT ll yield a slightly more complex interpretation: not only do 
intracomponent correlations have to be stronger than some in- 
tercomponentones but there are also conditions between inter- 
component coefficients. This may be more easily understood 
within the (3535, n22) state, where n is left as a variable, sim- 
ilarly to the SU(2) case discussed in Sec. IV A. Again, we 
have chosen this state purely for illustration reasons. 

Fig. |4]plots all filling factors and the relative signed eigen- 
value At,. As for the case K — 2, the graph is split into parts 
Lllandlll. 
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FIG. 4: Stability of the (3535, «22) wavefunctions. Filling factors v\ and V3 
(long dashed line), V2 and V4 (short dashed line), and the third eigenvalue /I3 
(solid line) of the (3535, «22) generalized Halperin wavefunction are plotted 
as a function of n. For n < 2 (part I) all quantities are positive and the corre- 
sponding state is stable. Part II corresponds to unphysical negative densities 
with a negative eigenvalue (-0.08). For n > 4, because of the eigenvalue /I3 
being still negative, the plasmas (l)-(3) phase separate from (2)-(4). 



• Part I (n < 2). One notices that type-(l/3) and (2/4) 
particles act identically therefore one can virtually treat 
the problem as for the case K - 2. Type-(l/3) and 
(2/4) plasmas are stable separately, with respect to the 
previous section on SU(2) wavefunctions. Moreover, 
type-(l/3) particles carry a V3 charge and are affected 
by type-(2/4) quasi-static impurities of nj V3 < V3 



charge. Similarly, type-(2/4) particles with charge V5 
interact with type-(l/3) plasma through an/ VB" < V5 
charge. Hence, one plasma suffers weaker repulsion 
from the other and type-(l/3) and (2/4) will mix in or- 
der to form a stable homogeneous state. 

• Part II (n = 3). The /I3 eigenvalue is now nega- 
tive (-0.08) as well as some filling factors. In the 
plasma picture, type-( 1/3) particles are equally repelled 
by their own species and type-(2/4) quasi-static impu- 
rities, each carrying a V3 charge. On the contrary, the 
type-(2/4) plasma still experiences more "favourable" 
repulsion from type-(l/3) particles. One cannot con- 
clude about the stability at this level and some care has 
to be taken of the inner composition of the type-(l/3) 
particles. Indeed, type-(l) particles cary a V3 charge 
and interacts with type-(3) quasi-static impurities via a 
2/ V3 charge which is less than any other charge for 
this plasma. Similarly, type-(3) will be less repelled by 
type-(l) particles than by any other particles. Hence, 
type-(l/3) plasma will tend to phase separate, which 
is contradictory to the phase mixing tendency of type- 
(2/4) particles. As in the K — 2 case, this yields un- 
physical negative densities for type-(l /3) particles. 

• Part III (4 < n). Above n = 4, all filling factors 
are positive but A3 becomes more and more negative. 
For n — 4, the same argument as above can be de- 
velopped: type-(l/3) plasma tend separate from type- 
(2/4), whereas type-(2/4) tend to mix with type-(l/3). 
Surprisingly, this is not related to a negative density, as 
in any other case previoulsy discussed. This may be 
due to the composite nature of type-(2/4) plasma. For 
n > 5, both plasma are more severely repelled by each 
other such that the phase separation is complete. 



C. Comparison with Exact-Diagonalization Studies 

The stability discussed so far is only related to the particu- 
lar form of the wavefunctions and it has somehow to be linked 
to the true ground-state of the quantum system with N inter- 
acting electrons. We therefore investigate the stability of gen- 
eralized Halperin wavefunctions within exact-diagonalization 
studies. The system is mapped onto a sphere in the center of 
which a magnetic monopole is fixed to ensure a magnetic field 
orthogonal to the surface. This magnetic monopole creates 25 
flux quanta threading the surface of the sphere. At a particular 
filling factor vj, the relation between the number of particles 
and that of flux quanta is 

2S = N/v T - 6, 

where 

Vt = YjMJ}, (17) 

u 
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and the shift 



hi 



(18) 



is due to the finite-size geometry and depends on the particular 
wavefunction considered. 

All calculations are performed within the lowest LL, with 
the help of Haldane's pseudopotentials, 36 VL which deter- 
mine the interaction between two electrons of type i and j 
with a relative angular momentum I. Halperin's wavefunc- 
tions (trii, rijj) represent the exact ground-state for a model in- 
teraction such that 



VI = 



vi 



1 for I < nii 

for / > nti 

1 for / < rtjj 
for / > riij 



(19) 



One of the simplest non-trivial cases occurs when all intra- 
(inter-)component correlations are the same, i.e. m, = m and 
«y = n. These states are fully unpolarized and the correspond- 
ing filling factor and shift are 



\m,n) 



K 



m + {K - \)n 



and St, 



(20) 



respectively. We have already shown in Sec. IV A and Sec. 
IV B that, within the plasma picture, m > n yields a stable and 
m < n an unstable state, whereas m = n represents a (stable) 
Laughlin state with SU(ZQ ferromagnetic order. This stability 
criterion may also be obtained directly from the interaction 
model corresponding to the (m, = m, n, 7 = n)-Halperin state: 
whenever this state is unstable, other zero-energy states with 
respect to the model interaction appear in the fully polarized 
sectors. These zero-energy states are quasihole excitations of 
the (mi = ny = m)-Laughlin state the model interaction of 
which matches that of the Haperin state in that sector. This is a 
direct consequence of v^"'" 1 < l/m for unstable states, as may 
be seen from Eq. (l2Qb with n > m. Thus any Zeeman-type 
perturbation or any extra pseudopotentials beyond the model 
interaction dramatically change the polarization, as suggested 
by the plasma picture. For more generic Halperin wave func- 
tions, similar conclusions can be drawn when one of the m, is 
lower than 1/vr- 

When the polarization is regarded as fixed, the phase sep- 
aration is clearly observed through the study of pair corre- 
lation function, which are discussed in Sec. V B. As men- 
tionned before, unstable systems will rather exhibit the insta- 
bility through unphysical correlation functions^ 

In addition to these general stability arguments, we inves- 
tigate via exact diagonalization the K — 2, (331) and (113) 
wavefunctions studied by Yoshioka et aJ.?° The first state is 
realized for the particular model VL = Vh = V?, = 1, all 
other potentials being zero, and 2S = 2N - 3. The second 
one is related to V° = = V 2 }[ = 1 and 2S = 2N - 1. For 
N - 6 electrons, and correspondingly 9 and 1 1 flux quanta, 
exact-diagonalization calculations yield an energy gap of 0.8 
for (331) as compared to 0.01 for (113). Hence, the (113) 



state has an energy gap which is almost two orders of magni- 
tude smaller than the characteristic energy, which is set to one. 
One may therefore expect that the (113) state is much less sta- 
ble than the (331) state, as indicated by the plasma analogy 
and the abovementioned argument. This is indeed the case as 
may be seen when other pseudopotentials are chosen non-zero 
in a perturbative manner. We choose, for this investigation, to 
vary continuously = Vh = V?, and V l E = VL from zero 
to one, in the case of the (331) state, and V\ = Vh - V], and 
V E = y|, for (113). Our exact-diagonalization results show 
that the unpolarized state described by the abovementioned 
wavefunctions (there is an equal number of spin f and X) is 
conserved only for the (331) case. Moreover, there is indeed 
an instability of the (113) state such that even a small perturba- 
tion in the pseudpotentials (Vj =0.1 for example) completely 
polarizes the state, in agreement with the general polarization 
argument given above. 

We now focus on the unpolarized sector although the sys- 
tem is no more in its ground-state in the (113) case. The 
unpolarized sector may be physically relevant for a system 
with constrained polarizations such as for a bilayer config- 
uration with equal densities in both layers. Fig. [5] presents 
phase diagrams with compressible and incompressible states 
for (331) and (113) varying states. The overlap between the 
exact ground-state and the Halperin wavefunctions is repre- 
sented by gradual shading. For the (331) case, it can be 
observed that there is a finite energy gap around the values 
V 3 = V E = of the exact model, and the overlap between 
the true ground state and the trial wavefunction remains large 
whenever the pseudopotentials V\ and V l E remain in this area. 
Further increase of the potentials leads to a gap collapse at rel- 
atively large values of the pseudopotentials. Hence, this state 
is stable even in the case of more realistic interactions beyond 
the model situation. 

In the case of the (113) wavefunction, the gap collapses 
rapidly even at small values of (« 0. 1). A rather subtle per- 
turbation would therefore completely change the system since 
a zero energy gap is incompatible with any FQHE. Outside 
the funnel-shaped area of compressible states, cf. Fig. [5jb), 
the overlap can be either quite small (« 0, at larger values of 
V l A and V|) or quite large (« 1, in the vicinity of the model 
situation). This indicates that the incompressible states are 
described by states with different symmetry. However, even 
though the system appears as stable up to relatively high V^, 
the gap remains small (0.01 to 0.1). Moreover we emphasize 
that, in the case of an unfixed polarization, the ground state is 
no longer in the unpolarized sector. 

Notice that the phase diagram in Fig. |5jb) is not generic 
for all unstable wavefunctions, but may be attributed to the 
pathologic model interaction of (1 13). General considerations 
on stability should only treat the polarization and pair correla- 
tion functions arguments. 

In the same manner as for the case K - 2, we compare the 
plasma picture with exact-diagonalization results for K = 4. 
For the special subset of matrices previously discussed, it has 
been checked numerically that unstable states are related to 
the existence of partially polarized zero-energy states with a 
lower number of flux quanta. As in the case K = 2, those 
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V. GROUND STATE PROPERTIES 



(b) 




V(intra)1 



FIG. 5: Phase Diagram for (a) 331 and (b) 113 ground-states - Exact- 
diagonalization were performed with N = 6 electrons and (a) 25 = 9 (b) 
25 = 11 flux quanta. An unpolarized state is assumed. Pseudopotentials 
V-ff = Vjj are noted V(intra) whereas Vfj, is V(inter). The overlap between 
Halperin's wavefuncions (exact ground-state of the unperturbed model) and 
the exact-ground state is plotted in gray scale. We also indicate the separa- 
tion between incompressible and compressible states (i.e. whether the ground 
state is in the L = sector or not) by a black line. 



Although the plasma picture is a powerful tool for the study 
of intrinsic properties of Laughlin and generalized Halperin 
wavefunctions, as shown in the previous section, it gives in 
itself no indication of the physical state chosen by the true in- 
teraction Hamiltonian. The Hamiltonian (O is obtained from 
a formal mapping of the wavefunctions to a corresponding 
plasma model, but it is not related to the original Hamilto- 
nian of interacting particles in the lowest or partially filled 
higher LL. Indeed, incompressible quantum liquids, which 
display the FQHE, are not found at all possible filling fac- 
tors for which one may write down a trial wavefunction. E.g. 
in the lowest LL, a Wigner crystal is energetically favor- 
able at v < 1/6.5, 38 in the first excited LL a succession of 
FQHE states and Wigner and bubble crystals 39 gives rise to 
a reentrant integral quantum Hall effect, 40 and in even higher 
LLs stripe phases 41,42 yield a highly anisotropic longitudinal 
transport. 43 Which of these competing phases is indeed cho- 
sen depends on the precise form of the true interaction of elec- 
trons in a fixed partially filled LL. 

It is however possible to express the energy of Laughlin's 
wavefunction in terms of the 3D Coulomb interaction poten- 
tial 'V(r) = e 1 /er and the pair correlation function 

g(r)cc j d 2 z 3 ....d 2 z N \V m (zi=0,z 2 = r,zi,...,z N )\ 2 , 

apart from a normalization constant, 

E = j dVV(r)[g(r)-l]. (21) 

The plasma picture may be of use here because the pair cor- 
relation functions for electrons and for plasmatic particles are 
the same, as may be seen from Eq. (f2]). The pair correlation 
function may be expanded as 28 

8 (z) = l-e^ 2 + f;im\ n e-^\ (22) 
/?= 1 * ' 

where the prime indicates a sum only over odd n (due to Fermi 
statistics), and the expansion parameters c„ vanish in the large- 
n limit. These expansion parameters are constrained in several 

|z|-<0 , 

respects. First, the short-range behavior g > \z\ implies 

that c„ = — 1 for n < m. Second, particular properties of 
the logarithmic potential in the plasma picture can be used to 
derive sum rules, which act as further constraints , 24 - 44 i 45 i 46 



states will be favored when any Zeeman-like perturbation is 
introduced. The phase transition predicted within the plasma 
picture is recovered. We checked this criterion for several par- 
ticular wavefunctions, such as (3333,233) and (3333,311), 
and it appears that wavefunctions with unstable correspond- 
ing plasma do polarize, partially or completely, in agreement 
with the classical stability discussed in Sec. IV B. For a given 
fixed polarization, phase separation is likely to be observed, 
as in the K — 2 case. 



A. Sum rules for SU(K) pair correlation functions 

For wavefunctions with SU(/T) symmetry, there are K(K + 
l)/2 pair correlation functions if all densities are well-defined, 
i.e. if Mk is invertible. The correlation function between type- 
(/) and type-(j) electrons is denoted by gij(z), and one may 
generalize the expression d22b to the case of ^-component 
wavefunctions, 

00 O/l |2\" 

S,( Z )= l_ e ^/ 2 + A ^ c ^e-^\ (23) 
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where the expansion coefficients c„ vanish for large n. Here, 
the prime indicates summation over odd n only for the intra- 
species functions g„. Indeeed, Fermi statistics is no more rel- 
evant when considering distinguishable electrons of type i and 
j + i. As for the Laughlin (K = 1) case, the short-range be- 



lz|-»0 , 

havior g u (z) > |zl " 



implies that c„ 



- 1 for n < n 



Further sum rules may be derived within the picture of K 
correlated plasmas introduced in the previous section. In or- 
der to derive those in the simplest manner, we decouple the 
different plasmas with the help of an orthogonal transforma- 
tion on the densities, 



(24) 



which diagonalizes the exponent matrix, Mk = P T DP, in 
terms of the orthogonal matrix P. The diagonalized Hamilto- 
nian thus reads 







>i(rr 




= P 











yv'v;(r)], 

ff, - - -"" 

I 



d 2 rjV P ;(r)-^ln|r-r'| P ;(r') 

2 



a An 



(25) 



where A/ is the i-th eigenvalue of M K and a t = £;[.P]y. The 
Hamiltonian < 7/[{p,'(r)}] is a sum of K independent Hamil- 
tonians 7f^[{p'(r)}], each of which corresponds to a single 
2DOCP. The correlation functions for these K plasmas must 
therefore obey the usual sum rules for 2DOCPs 24 ' 44 ' 46 



Mo = -1 
Mi 



M 2 = - 



M 3 



64 



-(1 - — ) 

(2ttM) 2 4 

(M-6X8-3M) 
(npAi? 



(26a) 
(26b) 

(26c) 

(26d) 



for the different moments 



M m = (p'y 



^rr^t^r)-!]. 



Here, primes indicate quantities in the diagonal basis. Eq. 
( I26at is due to the charge neutrality of the system, Eq. (I26bb 
reflects its perfect-screening property, and Eq. ( I26cl ) is a com- 
pressibility sum rule. The third moment [Eq. (I26dl )l has no 
apparent physical interpretation. Because the plasmas are de- 
coupled in the diagonal basis, there are no correlations be- 
tween different plasmas, i.e. g'^ir) = 1 for j + i. 

The pair correlation functions g,y(r) in the original basis 
may be obtained from the g^(r) with the help of the inverse 
orthogonal transformation. It is useful to start from the defini- 
tion of the structure factor in reciprocal space, which is related 
to the pair correlation function by Fourier transformation, 



5(k)-l=p j Sre^ [g(r)-l], 



for the simplest K — 1 case. It may also be expressed in terms 
of density operators, 



P f2S(k) = <p(k)p(-k)>-|<p(k)>| 2 , 



(27) 



where the quantities in brackets are averages with respect to 
the probability density function. In the case of K + 1, the 
structure factor has a matrix form, 



Q tffijStjOQ = <p,(k) Pj (-k)> - |<p,(k))| \(pj(k))\ 



and the associated pair correlation functions gy(r) may be 
obtained from those in the diagonal basis with the help of 
S y (r) = 6(r)6 u + (piPj) 1/z [gij(r) - 1], 



' P 2 [Sii(r)-l] ■ 


■ Pipdgiicir) - 


-ir 


= S(T) 




Pi ■ 


■ (T 


+ P T 


'p\ ■ 


■ (T 


P 


i 

i 

+P T 


■PKP\[gK\(r) - 1] ■ 


■ P 2 K lgKK(r) - 


i] , 






,o ■ 


' PK, 




,0 ■ 


' P'k> 




V 







(28} 



The sum rules ( I26al i to ( I26db for the diagonal basis thus im- 
mediatly yield those for the K correlated 2DOCP, due to Eq. 



For instance, the zeroth- and the first-moment rules are 



J d 2 rpi[ gij {Y) - 1] = -Sij 



(29a) 



J d 2 r^ Pi pj[ gi j(r)-l] = -^|.V/ ; (29b) 
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FIG. 6: Pair correlation functions related to (a) (331) for N = 10 particles 
and (b) (113) states for N = 8 particles, plotted as a function of distance (in 
units of magnetic length). Intra-component pairs (gn = #22) are plotted as 
a solid line and inter-component (gi2 = g2i) as dashed ones. In the (331) 
graph, all functions go to 1 at infinity, which is typical of uniform density. On 
the opposite, the (113) intra-component function vanishes at large distances 
thus endowing a particles aggregate. 



respectively. Eqs. j29aj and d29b| ) are generalizations of re- 
sults, which have previously been obtained for K = 2*&£L 
Second- and third-moment rules may be derived in the same 
manner. 



B. Exact-diagonalization results for pair correlation functions 

The pair correlation functions of the ground-state may also 
be obtained from exact-diagonalization results. We first focus 
on the SU(2) case by considering the same (331) and (113) 
wavefunctions discussed above and studied by Yoshioka et 
a/3 As mentioned in Sec. IV C, the first state is realized 
for the particular model VL = Vl, = V9, = 1, all other poten- 



TT ~ * II ~ ' U 

tials being zero, and 2S = 2N - 3. The second one is related 
to V T ° X = = V 2 }[ = 1 and 25 = 2N - 1. Fig. Eta) shows 
the ground-state pair correlation functions for the (331) state 
with N = 10 electrons and 25 = 17 flux quanta, and the re- 
sults for the (113) case (N = 8 and 2S - 15) are displayed 
in Fig. [6jb). The correlation function of a pair of electrons 



FIG. 7: Pair correlation functions related to (a) (3333,111) and (b) 
(3333,033) states for N = 8 particles, plotted as a function of distance 
(in units of magnetic length). Intra-component pairs (g„) are plotted as a 
solid line. In the (3333,111) graph, the dashed line corresponds to the inter- 
component pairs (gi£j, all being equal). In the (3333,033) graph, the dotted 
line shows the correlation function for electrons with opposite spin within 
any of each layer/valley. The dashed line corresponds to inter-component 
pairs between two different valleys or layers which are trivially constant in 
this non-correlated valley/layer example. 



(2, j) is the relative density of type- j electrons when a type- 
i electron is fixed at the origin, Pj(r) = pigij(r). With this 
definition in mind, one can infer that the (331) state is rather 
well-behaved. Correlation functions vanish near the origin, as 
a result of repulsive interactions and Fermi statistics. The den- 
sity peaks at finite distances indicate different layers of elec- 
trons (on average) with regular alternation of particles of type 
(1) and (2). The large-distance limit of type-(l) particles is 
related to a uniform density. Hence, the system discussed is 
composed of a uniform and homogeneous mix of electrons of 
both types. Unlike this first case, the (113) state displays a 
vanishing intra-component function (gn = #22) at large dis- 
tances. Electrons of both types thus tend to aggregate on a 
finite-size location. Moreover, the inter-component function 
(#12 = 821) is maximum only at infinity which implies that 
the different types of electrons tend to spatially separate. This 
corroborates the plasma picture of phase-separated particles. 
We now turn to the SU(4) case by considering the two sta- 
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ble generalized Halperin wave functions studied in Ref. l23l 
namely (3333,111) and (3333,033). The correlation functions 
are computed for N = 8 electrons and 2S — 9 flux quanta. 
While the partial filling factors of the (3333,111) state are 
fixed and all equal to 1/6, only the filling factor per layer is 
fixed for the (3333,033) ferromagnetic state. For a better com- 
parison, we have set the partial filling factors to be the same 
as those of (3333,1 1 1). Thus in both cases, there are only two 
electrons per species leading to prominent finite size effect 
(N = 12 is out of computational reach). The (3333,1 1 1) state 
is the straightforward SU(4) generalization of (331) and is 
therefore similar to its SU(2) counterpart (bearing in mind the 
small system size we are considering). The (3333,033) state 
consists of two independent spin insensitive v = 1/3 Laugh- 
lin states in each layer. Both intra-layer and intra-component 
correlation functions are therefore identical to their Laughlin 
counterpart (up to normalization factors). The inter-layer pairs 
are trivially constant in this non-correlated layer (3333,033) 
state. 

Furthermore, we have confirmed numerically the validity of 
the sum rules d29ab and d29bl ) for the SU(2) and SU(4) states 
discussed above. Notice that in the (3333, 033) state may be 
described alternatively by a SU(2) wavefunction, where the 
two components correspond to the two layer indices, regard- 
less of the spin orientation. The relevant sum rules are there- 
fore those of the SU(2) case, in which the exponent matrix Mk 
in Eq. d29b| ) is invertible. 

VI. CHARGED EXCITATIONS 

In this section, we study the charged excitations of SU(/T) 
Halperin wavefunctions with the help of the plasma analogy. 
For K = 1, quasi-hole excitations of the Laughlin wavefunc- 
tion may be written as 

N 

Y](zi-ZoW m ({z k }). (30) 
i=i 

They consist of adding a zero of density, or a magnetic flux, 
in the electron liquid at the position zo- 

Within the plasma picture, the extra Jastrow term adds a 
new potential term to the Hamiltonian ([3]), 

N 

tt^Ww-J^ta-zol. (31) 
i=i 

which thus describes a 2DOCP with a fixed impurity at zo 
and charge 1/ -\/m. Because of the plasma's perfect screening 
ability, 25 the particles are rearranged so that they screen the 
effect of the impurity in its vicinity. This requires l/m parti- 
cles in the plasma picture such that in the true electron liquid, 
the real charge e* of the excitation must be 

e = -, (32) 
m 

in units of the electron charge. Hence one electron can screen 
m excitations. 



We now investigate SU(ZQ Halperin wavefunctions for 
which there exist different types of excitation. The quasi-hole 
wavefunction 

N, 

n(z?-Zo)x^ ( .! ) m^ y (33) 

k,=\ 

creates here an excitation of (/)-type electrons at the position 
zo, i.e. adds a magnetic flux in the (0-th component of the 
electron liquid. The modified Hamiltonian in the plasma pic- 
ture contains a new potential of an impurity that only affects 
particles of type (j), 

^* = ^-5>|^- Z o| • 

Each of the K correlated 2DOCP exhibits perfect screening 
ability, i.e. the plasma of type (/) must screen the impurity 
totally, 

m,ef +Y J n ij ef = 1, (34) 

whereas plasmas of type (j), with j + i, screen a zero impurity 
^f+2>*ef = 0. (35) 

k 

Mi) 

Here, e ■ is the quasiparticle charge, in units of the electron 
charge, carried by electrons of type (j) in the electron liquid 
for excitations in the (i) component. Eqs. d34b and d35l) were 
previously derived for the K = 2 case, 6 and one may write 
them in a concise matrix form as 

Yj n ^T= S H ° ef=(M- K %. (36) 

The last equation is valid only if Mk is invertible. Indeed, 
if this is not the case, some component densities remain un- 
fixed, as described in the previous section, and one could only 
consider excitation of groups of particles with a definite den- 
sity. For instance, the (mmm) Halperin wavefunction has a 
non-invertible exponent matrix; the densities p\ and p2 may 
fluctuate although their sum remains fixed, vj — 1 /;«. Physi- 
cal excitation must therefore not distinguish between the two 
components. 

For the K — 4 case, (3333,1 1 1) wavefunction 23 exhibit four 
excitation types each of which carries a 1/6 charge, whereas 
for the (3333,033) wavefunction, only joined (l)-(3) and (2)- 
(4) excitations are to be considered, each of charge 1/3. 

Tab. U shows examples of charged excitations for SU(2) 
wavefunctions. The first example describes two independent 
Laughlin states and it is consistent that a type-(l) excitation 
should only affect type-(l) particles. The four next examples 
are related to the "(mmn) and (nnm)" problem. One should 
notice that there are inconsistencies for (1 13) and (223) states. 
Indeed, when an extra flux quantum is added to the type-(l) 
component, the number of flux quanta increases by one and 
the electron density remains the same. Therefore, the v\ filling 
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TABLE I: Charged excitations of Halperin's wavefunctions 

factor should decrease. However the total charge is conserved, 
so the sign of the quasi-hole charge should be the same as that 
of the electron, in order to compensate this "electronic" lack. 
This is not the case for (113) and (223) states which must be 
thus considered as unphysical, in addition to the conclusions 
drawn in the previous sections. The last example in Tab. |T]is 
simply a Laughlin wavefunction split into two arbitrary sets. 
There is only one common excitation, as for the usual U(l) 
case. 

Similarly, Tab. [TT] shows some charged excited states for 
SU(4) wavefunctions. All examples are associated with in- 
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e .(4) 


3333,111 
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1/6 


1/6 


1/6 


1/6 


3555,222 


2/5 


1/5 


1/15 


1/15 


1/15 


3535,222 


8/19 


3/19 


1/19 


3/19 


1/19 


5555,222 


4/11 


1/11 


1/11 


1/11 


1/11 



TABLE II: Charged excitations of SU(4) Generalized Halperin's wavefunc- 
tions 

vertible matrices Mk- No proof is given here, but we conjec- 
ture that unstable states yield some inconsistencies concerning 
the charged excitations, as for the K — 2 case. 

VII. CONCLUSIONS 

In conclusion, we have investigated the stability of Halperin 
wavefunctions for /^-component quantum Hall systems, with 
a particular emphasis on the cases K — 2 and 4. The asso- 
ciated SU(2) and SU(4) internal symmetries happen to be the 
physically most relevant if one considers, e.g., bilayer quan- 
tum Hall systems and graphene in a strong magnetic field. The 
K — 4 case occurs when the Zeeman effect is relatively small 
with respect to the leading interaction energy scales. In or- 
der to derive the stability criteria, we have generalized, in a 
systematic manner, Laughlin's plasma analogy to multicom- 
ponent systems. The validity of the criteria is corroborated 
with the help of exact-diagonalization studies. 

As for the conventional one-component quantum Hall sys- 
tem, the quantum-classical analogy yields a compelling phys- 



ical interpretation of the trial wavefunctions, in terms of K 
correlated 2DOCP. Besides the stability of the trial wavefunc- 
tions, it also allows one to understand relevant ground-state 
properties, such as the associated pair-correlation functions, 
and fractionally charged quasiparticle excitations. 

Whether the discussed trial wavefunctions correctly de- 
scribe the true ground state in physically relevant multicom- 
ponent systems, such as bilayer quantum Hall systems or 
graphene, depends on the precise form of the interaction po- 
tential. The plasma analogy with its rather artificial inter- 
action may not give insight here, and variational or exact- 
diagonalization studies need to be performed to determine the 
correct ground state for a physical interaction potential. It 
has indeed been shown that in the case of Coulomb interac- 
tion, a possible FQHE state at v = 2/3 is not described by 
a generalized SU(4) Halperin wavefunction. 22 ' 23 More com- 
plicated trial wavefunctions, such as composite-fermion 2 ^ or 
even more exotic states, may describe FQHE states at this 
and possibly other filling factors in a more appropriate man- 
ner. However, in the U(l) one-component quantum Hall sys- 
tem, the inevitable starting point in the understanding of the 
FQHE is Laughlin's wavefunction; 2 other wavefunctions may 
be viewed as sophisticated generalizations of it. In the same 
manner, the study of SU(.K) Halperin wavefunctions and the 
plasma analogy yield important physical insight into multi- 
component quantum Hall systems, and one may conjecture 
that they play a similar basic role for possible generalizations 
as Laughlin's in the U(l) case. 

Furthermore, it has been shown in the SU(2) case, that 
not all possible, though stable from our analysis, Halperin 
wavefunctions are valid candidates from a symmetry point 
of view. Indeed, most of the (m,m,n) wavefunctions are 
not eigenstates of the total spin operator, the Casimir oper- 
ator of SU(2), as they should for spin-independent interac- 
tion Hamiltonians. 47 The (331) wavefunction discussed here, 
is, e.g., not an eigenstate of the total spin. However, this 
problem may be cured by attaching the permanent of the ma- 
trix (z^ - z^Y 1 to the (331) wavefunction. 47 ' 48 The situa- 

' J 

tion is more complicated in the SU(4) case, where there are 
more Casimir (spin-pseudospin) operators, and the wavefunc- 
tions should be eigenstates of these operators. More detailed 
theoretical investigations are required to settle the question 
whether some SU(4) wavefunctions may be corrected in a 
similar manner. 
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